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ABSTRACT 

We present multi-frequency VLA observations of the radio galaxy PKS 1246-410 at 
the center of the Centaurus Cluster, and compare these results to recent Chandra 
observations. The unusual radio morphology of PKS 1246-410 probably results from 
interactions with the hot, X-ray emitting gas. This gas, along with cluster magnetic 
fields, also produces substantial Faraday Rotation Measures towards PKS 1246-410. 
We discuss these observations in the context of a sample of 14 luminous X-ray clusters 
with embedded radio galaxies and possible cooling flows. A correlation is found be- 
tween the cooling flow rate and the maximum Faraday Rotation Measures. Magnetic 
fields of strength 10-40 fiG are found to be common to the centers of clusters with 
strong cooling flows, and somewhat lower field strengths of 2-10 [iG are found in the 
non cooling-flow clusters. 

Key words: galaxies: active - galaxies: nuclei - radio continuum: galaxies; galaxies 
individual - PKS 1246-410 



^ ■ 

Or 
i 

o 

to 

in 



1 INTRODUCTION 

Radio sources embedded within rich clusters of galaxies can 
serve as useful probes of the intracluster medium (ICM). 
The thermal ICM pressure may confine the radio source and 
polarized radiation from embedded radio sources may also 
be rotated by the Faraday effect if magnetic fields are present 
in the ICM. 

The X-ray emission in many clusters is strongly peaked 
at the center, indicating high densities and cooling times of 
the hot ICM in the inner ~100 kpc which are much less 
than the Hubble time. To maintain hydrostatic equilibrium, 
either another source of pressure or an inward flow is re- 
quired (for a review see Fabian 1994). The cooling flow phe- 
nomenon in clusters is quite common with the short cooling 
time condition occurring in 70-90% of an X-ray flux-limited 
sample of clusters (Edge, Stewart & Fabian 1992; Peres et 
al 1998). 

While most extragalactic radio sources exhibit Faraday 
rotation measures (RMs) on the order of 10s of rad m -2 
due to the interstellar medium of our galaxy, there are a 
small number of sources displaying RMs in excess of ~1000 
rad m~ 2 . Studies of these high RM sources (e.g. Cygnus A: 
Dreher, Carilli & Perley 1987; A1795: Ge & Owen 1993; 
Hydra A: Taylor & Perley 1993) indicate that these RMs 
are most likely to be induced by a cluster magnetic field. 
The high RM sources are typically surrounded by luminous, 
X-ray emitting gaseous halos, and in most cases a cooling 
flow (Taylor, Barton & Ge 1994). The diversity of the radio 
structures and powers, and the commonality of the X-ray 
parameters suggest that the high RMs are the result of ex- 



ternal Faraday rotation occurring on large scales, and are 
not due to a local interaction between the radio source and 
its environment. For external Faraday rotation, the RMs are 
related to the density, n e , and magnetic field along the line- 
of-sight, Bu, through the cluster according to 



RM = 812 



c _B|idi radians m 



(1) 



where Bn is measured in /iGauss, n c in cm and dl in 
kpc. The RM distribution, along with X-ray observations 
(used to estimate n c ), can then be used to understand the 
magnetic field structure along the line-of-sight (see review 
by Carilli & Taylor 2002). 

Here we present VLA (Very Large Array )^ observations 
of the heretofore little known radio galaxy PKS 1246-410 
at the centre of the Centaurus cluster. This radio galaxy 
should not to be confused with the more famous radio galaxy 
Centaurus A. In §2 we describe the radio observations. In §3 
we summarize the X-ray observations, and in §4 we present 
our analysis of both the radio and X-ray observations. In 
§5 we present our current knowledge of the sample of 14 
clusters with embedded radio galaxies originally selected by 
Taylor et al. (1994). 

We assume Ho = 50 km s _1 Mpc -1 and ft = 1 through- 
out. 



* The National Radio Astronomy Observatory is a facility of the 
National Science Foundation operated under a cooperative agree- 
ment by Associated Universities, Inc. 
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2 RADIO OBSERVATIONS AND DATA 
REDUCTION 

Observations of PKS 1246-410 were made at 1.465, 1.565, 
4.635, 4.885, 8.115, and 8.485 GHz at the VLA. The de- 
tails of the observations are provided in Table 1. The source 
3C 286 was used as the primary flux density calibrator, 
and as an absolute reference for the electric vector polar- 
ization angle (EVPA). Phase calibration was derived from 
the nearby compact source J1316-3338 with a cycle time 
between calibrators of about 12 minutes. Instrumental cali- 
bration of the polarization leakage terms was obtained using 
OQ 208 which was observed over a wide range in parallactic 
angle. The data were reduced in AIPS (Astronomical Image 
Processing System) , following the standard procedures. The 
AIPS task IMAGR was used with a suitable taper to make 
/, Q and U images at each of the 6 frequencies observed 
at the same resolution. Polarized intensity, P, images and 
polarization angle, X-i images were derived from the Q and 
U images. The Faraday Rotation Measure (RM) image was 
made from the combination of the 4 higher frequency \ lra ~ 
ages for PKS 1246-410. Pixels in the RM image were flagged 
when the error in any of the \ images exceeded 15 degrees. 
Stokes / images using multiple frequencies within an ob- 
serving band were also produced. In addition to the added 
sensitivity, the image produced also benefits from improved 
(u, v) coverage. No correction for the spectral dirty beam 
has been performed, but the sidelobes expected from this 
effect should be below the thermal noise floor. 



3 X-RAY OBSERVATIONS 

Allen et al. (2001a; see also Ikebe et al. 1999) present ASCA 
observations of the Centaurus cluster. They find that the 
spectrum of the central 100 kpc can be well-described by 
either a two-temperature or cooling flow model with a mass 
deposition rate of 30-50 Mq yr _1 . 

Sanders & Fabian (2001) present data from Chandra 
observations at much higher spatial resolution and report 
the discovery of a knotty, plume-like feature that begins 60" 
(17 kpc) to the northeast of the nucleus, and curves to the 
south before terminating within 1" (0.28 kpc) of the radio 
core (see §4.2). Sanders & Fabian find that this plume is 
0.25 keV cooler than the surrounding thermal gas and report 
on the presence of a "cold front" in the Centaurus cluster 
where the temperature sharply increases near the end of the 
plume, to the north-east of the core. One way to produce 
such features is by suppression of thermal conduction by 
magnetic fields (see Carilli & Taylor 2002 and references 
therein). 

The electron density profile for the central 70 kpc (~ 
4.2') radius region, determined from a deprojection analysis, 
is shown in Fig. 1. The best-fitting /?— model for the density 
has parameter values no = 0.09 ± 0.004 cm -3 , (3 = 0.62 ± 
0.02, and a core radius, r c = 5.66 ± 0.04 kpc. Note that the 
core size is similar to that of the radio source, which has 
made 'holes' in the X-ray emission at radii of 1-5 kpc. 



4 RESULTS 



4.1 Radio Properties of PKS 1246-410 

The radio source PKS 1246-410 is associated with 
NGC 4696, an elliptical galaxy at the center of the Cen- 
taurus cluster (also known as Abell 3526). The redshift is 
0.0099, so 1" corresponds to 0.28 kpc. The radio power 
of PKS 1246-410 is 1.5 xlO 24 W Hz" 1 at 1.6 GHz, iden- 
tifying it as a low power radio galaxy for which we would 
expect an FR 1 morphology (edge dimmed "woolly" radio 
tails). In fact the radio morphology at 1.6 GHz, shown in 
Fig. 2, appears quite different from most low power radio 
galaxies. The usual well defined core and thin jets feeding 
into larger lobes are missing. Instead there is a bright point 
source which we identify as the core based on its central lo- 
cation and slightly flatter spectral index of a = —0.5, where 
S„ oc v a (Fig. 3). The position of the core and a summary 
of source properties is given in Table 2. The core is sur- 
rounded by "lobe" emission with a predominantly east-west 
orientation. At their extremes both lobes appear (at least in 
projection) bent to the south and back in towards the core, 
although there is also a weak extension to the northeast. The 
emission in the extremities is quite steep (a ~ —2). Over- 
all the impression is of a rather distorted radio morphology, 
perhaps due to confinement by the hot X-ray emitting gas. 
Similar distortions have been seen for other radio galaxies 
embedded in dense X-ray environments with the most ex- 
treme case being the amorphous radio galaxy PKS 0745-191 
(Taylor et al. 1994). The swept-back appearance of the outer 
parts of the lobes may be due to relative motion of the galaxy 
and ICM (see Sanders & Fabian 2001 for more discussion on 
this point). 

Fortunately, and unlike PKS 0745-191, the bright radio 
lobes of PKS 1246-410 are moderately well polarized (5 - 
40%) , so it is possible to derive the RM distribution across 
nearly the entire radio source. This has been done using 
four frequencies spread across the 4 and 6 cm bands, and 
is shown in Fig. 4. The RMs appear patchy and are not 
correlated with any features in total intensity. The RMs in 
PKS 1246-410 are typically in the range -500 to +500 rad 
m -2 and reach a maximum of 1800 rad m~ 2 . A histogram 
of the RM distribution is shown in Fig. 5. After correcting 
for Faraday Rotation, the projected intrinsic magnetic field 
orientation of the synchrotron emission generally follows the 
edges of the source (Fig. 6) . This is typical for extended radio 
sources. 



4.2 X-ray/Radio Comparisons 

In Figs. 7 and 8 we show overlays of the X-ray and radio 
emission. This cluster shows a dramatic plume in the X- 
rays that reaches in to the cluster center and wraps around 
the radio galaxy. At the highest spatial resolution the radio 
emission seems to emerge from the core along a nearly N- 
S axis (Fig. 9). The rapid change to a more nearly E-W 
orientation of the radio emission is most likely the result of 
interactions with the thermal gas. Further from the nucleus 
at ~3 kpc the radio lobes bend to the south. The eastern lobe 
appears to fill in the region bounded by the plume, while the 
western lobe also spills around it. In general the impression 
is that the radio source has been contained by the dense 
gas in the plume. Similar examples of radio emission filling 
X-ray cavities and avoiding overdensities have been found 
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for 3C 84 at the center of the Perseus cluster (Bohringer et 
al. 1993, Fabian et al. 2000), in PKS 1508+059 in A2029 
(Taylor et al. 1994) and in Hydra A in A780 (McNamara et 
al. 2000). 



4.3 Magnetic Field Estimates 

The gas density profile derived from ROSAT HRI observa- 
tions can be parameterized with a /3-model fit as described 
in §3. This gives an electron density at radius 10 kpc of 0.04 
cm~ 3 . Considering the patch of —1500 rad m~ 2 located 3 
kpc west of the core, the minimum constant magnetic field, 
B\\ , that could produce the observed RM over a path length 
of 10 kpc has a strength of 8 /iG. Several changes in sign 
are in evidence on both sides of the core, however, indicat- 
ing that the magnetic field direction is rather chaotic. The 
simplest explanation for the RM distribution is a tangled 
magnetic field with a coherence length of ~1 kpc. For a 
density distribution that follows a /3-profile, Felten (1996) 
has derived the following relation for the RM dispersion: 



CRM = 



KB n rl /2 l 1/2 / r(3/3-0.5) 
(l + r 2 /r2)(6/3-i)/4 y r(3/3) 



(2) 



where no, r c , and (3 are derived in §3.2, I is the cell size, r is 
the distance of the radio source from the cluster center, Y is 
the Gamma function, and K is a factor which depends on the 
location of the radio source along the line-of-sight through 
the cluster: K = 624 if the source is beyond the cluster, 
and K = 441 if the source is halfway through the cluster. 
Note that Eq. 2 assumes that the magnetic field strength, 
B, is related to the component along the line of sight, (By), 
by B = %/3-Bii . Given the clear interactions between the X- 
ray and radio plasma, we can assume that PKS 1246-410 is 
close to the true center of the cluster, and not fortuitously 
projected on the center. With the above measurements or 
assumed values, we use the estimated RM dispersion of 660 
rad m -2 from Fig. 5 to estimate the cluster magnetic field 
strength at 11 uG. In §5.2 we compare this value to those 
for a larger sample of clusters. 

Fig. 8 shows how the X-ray emission varies across the 
radio lobes. It is interesting to compare it with Fig. 4 which 
shows the rotation measure map. The parts with positive 
RM (blue in Fig. 4) near the centre have the highest X-ray 
surface brightness, <5x (brown in Fig. 8); the parts with neg- 
ative RM (yellow/orange in Fig. 4) have a lower Sx (green 
in Fig 8) and the parts with little RM (green in Fig. 4) 
have the lowest Sx (blue in Fig. 8). There is thus tentative 
evidence that the regions of high \RM\ have high Sx and 
at least some of the excess X-ray emission lies in front of 
the radio source. Since \RM\ oc nBi and the X-ray surface 
brightness Sx oc n 2 T 1 ^ 2 ^ (the dependence on temperature 
T is actually complicated by line emission so the propor- 
tionality is approximate), if the magnetic field B oc VnT, 
i.e. the magnetic pressure is in equipartition with (or a fixed 
fraction of) the thermal pressure, then \RM\ oc SSxn^ 1 ^' 2 , 
where SSx is that part of Sx which is in front of the radio 
source. Qualitatively this explains what is seen, but the de- 
tailed relationship must depend upon the field configuration, 
which must be complex. 



£ 02 




5x10 



0.01 0.02 
Radius (Mpc) 



Figure 1. Dcprojected density profile determined from the Chan- 
dra observations. The solid line shows the best-fitting /3— model. 
The inner few kpc are affected by the 'hole' in X-ray emission 
associated with the radio lobes. 
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Figure 2. VLA image of PKS 1246-410 at 1.565 GHz. The 
greyscale has been chosen to show the faint extended emission 
in the range to 5 mjy/beam. Contours start at 0.3 mjy/bcam 
and increase by factors of 2. The peak flux density in the map is 
358 mjy/beam. The restoring beam has dimensions 4.4 X 2.1" 
in position angle degrees. 
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Table 3. Radio Sources in Cooling Flow Clusters 
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a The mass flow rate for a critical cooling time of 1.3 X 10 10 from Allen & Fabian (1997), White, Jones &; Forman (1997) and Peres 
et ai. (1998). 

b Galaxy type (E: elliptical; pec: peculiar, or cD: central dominant elliptical) 
c The largest angular size of the radio source in kpc. 

References - (1) Dreher et al. 1987; (2) Owen et al. 1990; (3) This work; (4) Burns et al. 1992; (5) Taylor & Perley 1993; (6) Ge & Owen 
1993; (7) Feretti et al. 1999; (8) Taylor, Barton, & Ge 1994; (9) Taylor, in prep.; (10) Taylor et al. 2001 
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Figure 3. A spectral index map for PKS 1246-410 made by com- 
bining images at 4.76 and 8.30 GHz. Contours represent the 4.76 
GHz total intensity starting at 0.25 mjy/beam and increasing by 
factors of 2. The peak in the image is 106 mjy/beam. The restor- 
ing beam is plotted in the lower right corner and has dimensions 
2.1 X 1.2 in p.a. 19°. 
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Figure 4. The rotation measure (RM) map for PKS 1246-410 at 
the center of the Centaurus cluster. Contours are the same as in 
Fig. 3. 
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Figure 5. Rotation measure distribution for PKS 1246-410. The 
mean value is —356 rad m -2 with a dispersion of 662 rad m — 2 . 
The width of each bin is 8.3 rad m — 2 , which is about the same 
as the error in the RM determination. 



Figure 7. The Centaurus Cluster in the X-ray (false colour from 
the Chandra image) and radio (8.4 GHz contours). The resolution 
is similar in both images. 
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Figure 6. The Projected magnetic field vectors at 4.9 GHz (vec- 
tors corrected for Faraday rotation) overlaid on total intensity 
contours for PKS 1246-410. Contours are the same as in Fig. 3. 
The length of the vectors is proportional to polarized intensity 
(1" = 0.31 mjy/beam). 




Figure 8. The inner part of the Centaurus Cluster in the X- 
ray (ct arcsec -1 31.5 ks" 1 (0.5-5 keV); false colour from the 
adaptively-smoothed Chandra image) and radio (5 GHz con- 
tours). Contours are the same as in Fig. 3. 
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Figure 9. VLA image of PKS 1246-410 at 8.30 GHz made from 
the BnA configuration data. Contours start at 0.35 mjy/bcam 
and increase by factors of \/2. The peak flux density in the map 
is 63 mjy/beam. The restoring beam has dimensions 1.05 X 0.58" 
and a position angle of 21 degrees. 
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Figure 10. The maximum absolute RM plotted as a function 
of X for the 13 clusters in Tabic 3 with measured RMs. We also 
include the well-studied, but higher redshift 3C 295 cluster on this 
plot. 

5 A SAMPLE OF RADIO SOURCES 
EMBEDDED IN CLUSTERS 

5.1 Summary of results 

Examination of a flux-limited, all-sky X-ray sample of galax- 
ies put forth by Edge, Stewart, & Fabian (1992) reveals that 
clusters with non-zero mass flow rates frequently display sig- 
nificant RMs. Taylor, Barton & Ge (1994) formed a sample 



Table 1. Observational Parameters 



Source 


Date 


Frequency 
(GHz) 


Bandwidth" 
(MHz) 


Config. 


Duration 
(hours) 


PKS 1246 


Aprl99S 


1465/1565 


25 


A 


0.24 


-410 


Junl998 


1465/1565 


25 


BnA 


0.99 




Aprl998 


4635/4885 


50 


A 


0.09 




Junl998 


4635/4885 


50 


BnA 


1.23 




Novl998 


4635/4885 


50 


CnB 


0.61 




Junl998 


8115/8485 


50 


BnA 


0.89 




Novl998 


8115/8485 


50 


CnB 


0.62 



Table 2. Source Properties 



Property 


PKS 1246-410 


core RA (J2000) 


12 h 48 m 49f273 


core Dec. (J2000) 


-41° 18'39!Tl 


velocity 


2958 ± 15 km s" 1 


projected distance 


Mpc 


angular size 


122" 


physical size 


34 kpc 


flux density (5 GHz) 


1330 mjy 


power (5 GHz) 


5.7 X 10 23 W Hz" 1 


<RM> 


-356 rad m~ 2 


CRM 


662 rad m~ 2 


|RM max | 


1800 rad m" 2 



of 14 cooling flow clusters each containing at least one em- 
bedded radio galaxy stronger than 100 mjy at 5 GHz. All 
of these sources have now been extensively studied with the 
VLA, and in most cases considerably better X-ray observa- 
tions are also available. Here we present an analysis of the 
combined results for this sample. 

Of the 14 clusters in the sample, 10 display an excess 
of 800 rad m~ 2 , two have RMs between 300 and 700 rad 
m~ 2 (one is A119 discussed by Feretti et al. 1999; and the 
other is the 3C 129 cluster), and two could not be measured 
(PKS0745-191 and the Perseus cluster) due to a lack of 
polarized flux. Nominal mass-flow rates, X-ray luminosities 
(from 2-10 keV), RMs, and other relevant parameters for 
these 14 clusters are listed in Table 3. Recent work has shown 
that the nominal X-ray mass flow rates, calculated under the 
assumption that the cooling flows extend out to the cooling 
radii in the clusters (where t coo \ — 13 Gyr) are likely to 
overestimate the true mass deposition rates from the cooling 
flows by a factor of 3 or more. This is required by the fact 
that the cooling flows are much younger than 13 Gyr (Allen 
& Fabian 1997; Allen et al. 2001a) and, importantly, by the 
apparent lack of any gas cooling below 1-2 keV in spectra 
taken with XMM-Newton (Peterson et al 2001; Tamura et 
al 2001). For display purposes, we shall continue to use the 
nominal X-ray derived values here, simply as a characteristic 
measure of the amount of cooler gas in the cluster core. 
However, to signify that they are not complete mass cooling 
rates (although they compare well with the rate at which 
gas cools from the cluster virial temperature to about one 
third that value), we denote them as X instead of M. 

The RM value given is chosen from the region of largest 
RM which is coherent across a beamwidth or more with good 
polarized signal-to-noise. We have also plotted in Fig. 10 the 
RM as a function of cooling rate, X. These results strongly 
suggest a correlation between the cooling rate and rotation 
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measure. Even more striking is the fact that all of the known 
high RM sources (RM > 700 rad m~ 2 ; see Table 1 in Tay- 
lor, Inoue, & Tabara 1992) at low redshifts (z < 0.4) are 
included in this X-ray selected sample. The high RM ra- 
dio galaxy 3C 295 (Perley & Taylor 1991) at a redshift of 
0.461 is also associated with a luminous X-ray cooling flow 
(Henry & Henriksen 1986; Allen et al. 2001b), but the pre- 
Chandra X-ray data for the remaining seven high RM, high 
z, sources are not good enough to distinguish between ex- 
tended thermal emission and that from an active nuclear 
source. Beyond redshifts of ~0.4, low power radio sources 
like PKS 1508+059, and PKS 1246-410 are much too weak 
to allow for a determination of their RM with the current 
generation of radio telescopes. 



5.2 Cluster Magnetic Field Strengths and 
Topologies 

Our analysis of the RM dispersion in PKS 1246-410 leads to 
an estimated cluster magnetic field strength of 11 /iG in the 
Centaurus cluster. This is about a factor of two larger than 
the field strengths of 6 fiG estimated for the two non-cooling 
flow clusters in our sample (3C 129 and A119) by the same 
method. At the other end of the range, applying this method 
to Hydra A leads to an estimated field strength of 35 fiG. 

The scale size of the RM fluctuations in Centaurus is 
quite small at ~1 kpc. In cluster radio sources with both 
lower and higher RMs (e.g. A119, 3C129, Hydra A) typical 
scale lengths are ~3 kpc. 

The small physical size of PKS 1246-410 only allows 
us to probe the magnetic fields along the line-of-sight to 
the center of the cluster. Strictly speaking the region across 
which the RMs can be measured is only a projected dis- 
tance of ~10 kpc, so the RMs could be explained in terms 
of the host galaxy. However, in other clusters such as A119 
(Feretti et al. 1999) and A514 (Govoni et al. 2001), RM 
studies have been performed towards several radio galaxies 
either in or behind the cluster and the magnetic fields have 
been found to prevail to large radii, ~500 kpc. In a study 
of background radio sources seen through clusters, Clarke, 
Kronberg & Bohringer (2001) argue convincingly in favor 
of cluster-wide magnetic fields distributed throughout the 
central ~1 Mpc. 



6 CONCLUSIONS 

The radio galaxy PKS 1246-410 appears to be strongly influ- 
enced by the dense X-ray emitting cluster gas that surrounds 
it, especially the plume discovered by Sanders & Fabian 
(2001). The radio - X-ray interaction produces an unusu- 
ally steep spectrum, small source with bent lobes. Combined 
with a cluster magnetic field, the hot gas also produces Fara- 
day RMs reaching 1800 rad m~ 2 . With some assumptions 
about the field topology and distribution, a cluster magnetic 
field of ~11 fiG is derived. This field strength is intermedi- 
ate between that of clusters in which there is no cooling-flow 
present, and those estimated to have a large cooling- flow. 

In the near future the Expanded VLA will provide an 
order of magnitude increase in sensitivity over the current 
VLA. This should allow RM studies towards more than a 



tenfold increase in the number of sources within a given clus- 
ter. By analyzing the RM distribution at many lines of sight 
through a cluster it should be possible to determine how the 
magnetic field changes with cluster radius and density. 
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